Rotational motion of magnon and thermal Hall effect 
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Due to the Berry curvature in momentum space, the magnon wavepacket undergoes two types of 
orbital motions in analogy with the electron system: the self-rotation motion and a motion along 
the boundary of the sample (edge current). The magnon edge current causes the thermal Hall effect, 
and these orbital motions give corrections to the thermal transport coefficients. We also apply our 
theory to the magnetostatic spin wave in a thin-film ferromagnet, and derive expression for the 
Berry curvature. 

PACS numbers: 85.75.-d, 66.70.-f, 75.30.-m, 75.47.-m 



I. INTRODUCTION 

Recently in the spintronics field, the magnon (spin 
wavei) transport in an insulating magnet attracts much 
attention as a candidate of a carrier of the spin informa- 
tion with good coherence and without dissipation of the 
Joule heating. In particular, magnon can propagate over 
centimeter distance in some magnets^, e.g. yttrium- iron- 
garnet (YIG), and this is long enough compared with the 
spins in metals and doped semiconductors. The magnon 
current can be experimentally generated by the spin Hall 
effect^, and its motion can be observed by the time- and 
space-resolved measurement methods'*. Besides, a pre- 
cise control of the spin information is necessary for the 
application in spintronics devices. 

The thermal Hall effect (Righi-Leduc effect) of the 
magnon, which is useful to control the magnon trans- 
port, is predicted theoretically by Katsura et al.^ and 
observed experimentally by Onose et air . Katsura et al. 
considered a ferromagnet with the Kagome lattice struc- 
ture, and calculated a thermal Hall conductivity by the 
Kubo formula; Onose et al. measured the thermal Hall 
conductivity using an insulating ferromagnet LU2V2O7 
which has a pyrochlore structure with the Dzyaloshinskii- 
Moriya (DM) interaction. On the other hand, we found 
in our recent studyi that there are correction terms to the 
thermal Hall conductivity in the linear response theory, 
and showed that the thermal Hall effect of the magnon 
arise from the edge current of the magnon in the semi- 
classical picture. Our theory is applicable not only to the 
quantum mechanical spin wave, e.g. in LU2V2O7, but 
also to the classical magnetostatic spin wave, where the 
wave length is long and exchange coupling is negligible, 
e.g. in YIG film. 

In the present paper, we develop the transport theory 
of magnons with detailed calculations. Some parts of the 
theory has been published in Ref. 0- There are two ap- 
proaches; the semiclassical theory and the linear response 
theory. From the semiclassical equation of motion, the 
magnon edge current are described by the Berry curva- 
ture in momentum space and does not depend on the 
details of the system such as the shape of the boundary 
of the sample. From this magnon edge current, we obtain 
the magnon current and energy current density under a 



spatial variation of the temperature or the chemical po- 
tential, resulting in the thermal Hall effect of the magnon. 
Since our result of the thermal Hall conductivity does 
not agree with the previous works in Refs. 5,6, we refor- 
mulated the linear response theory in analogy with the 
electron system^, by noting that the temperature gradi- 
ent is not a dynamical force but a statistical force. It 
is identified that the difference from the previous work 
arises from orbital motions of the magnon, and that the 
magnon rotates around itself besides the magnon edge 
current. 

We apply our theory to LU2V2O7, and calculate the 
orbital angular momenta of the rotational motions of the 
magnon. For another application, the expression of the 
Berry curvature for the magnetostatic forward volume 
wave in YIG is derived. In this case the Berry curva- 
ture of the highest energy band enhances and that of the 
other bands converges to at A: = 0. Besides the Berry 
curvature becomes larger as the magnetic field becomes 
small. 

This paper is organized as follows. We present the 
semiclassical theory for the magnon and consider the 
thermal Hall effect of the magnon in Section [Hi The 
linear response theory with a temperature gradient and 
the orbital motions of the magnon are discussed in Sec- 
tion mil Section IIVI and Section |V] are devoted to ap- 
plications of our theory to LU2V2O7 and YIG, respec- 
tively. We conclude with a summary in Section IVIl and 
a brief review of the linear response theory for the elec- 
tron system and some useful equations are presented in 
Appendix |^ 

Throughout this paper we consider two-dimensional in- 
sulating magnetic systems for simplicity, and assume that 
magnons do not interact with each other. Generalization 
to three-dimensional magnets is straightforward. 



II. SEMICLASSICAL THEORY 

Our approach is based on the semiclassical theory, 
in analogy with the electron systent^i^. We consider a 
magnon wavepacket which is well localized around the 
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center (rc,fcc) in the phase space: 
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where \4>nk) is the Bloch wave function in the nth magnon 
band, an{k,t) satisfies 



J dk \anik,t)f = 1, 
J dk\an{k,t)\'^ k = kc. 



and \Wn) satisfies 



{Wn\r\Wn) 



(2) 
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(4) 



Hereafter we omit the index c for brevity. The dynam- 
ics of the wavepacket is described by the semiclassical 
equation of motion, which includes the topological Berry 
phase term: 



. _ 1 dSnk _ 
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Here n is the band index, Snk is the energy of the magnon 
in the rith band, n„(fc) is the Berry curvature in momen- 
tum space: 
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with \un{k)) being the periodic part of Bloch waves in 
the nth band defined as (/)„fc(r) — Un{k,r)e'^^''^ . U{r) is 
a confining potential which exists only near the boundary 
of the sample. This potential U{r) forbids the magnon 
wavepacket going outside of the sample, and its gradient 
exerts a force on magnons. Such approach of the confin- 
ing potential is successful in describing the edge picture 
of the quantum Hall effect in electron systemsiS. Thus 
we have similarly introduced the confining potential for 
magnons. Strictly speaking, for the validity of Eqs. ([5]) 
and dni, the spatial variation of U{r) should be much 
slower, compared with the size of the wavepacket. Nev- 
ertheless, as we can see from the quantum Hall effect as 
an example, many of the results for the slowly varying 
U{r) are expected to carry over to the case of rapidly 
changing U{r) as well. 

Near the edge of the sample, there exists an edge cur- 
rent of magnons due to the anomalous velocity term 

k X 0„(fc) = VU{r)/h X 0„(fc) in Eq. ©. For ex- 
ample, the magnon edge current for the edge along the y 
direction is expressed as 

ly^ f dx:i^p(£„fc + C/(r))[Vf/(r)/Sx a.(fc)]j,, 

•^^ n,k 

1 r°° 

= -wT. depie)n^Ak), (8) 
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FIG. 1: Coordinate of the ferromagnet, used for the calcula- 
tion of the edge current. U{r) is a confining potential. 



where x — a and x = b are chosen well inside and outside 
of the sample so that U{a) = and U{b) = oo, is the 
area of the sample, p{s) is the Bose distribution function 
p(e) = (e/^C^-M) _ p = l/knT, is the Boltz- 

mann constant, /i is the chemical potential, and T is the 
temperature. Henceforth the magnon current means the 
current of the magnon number. We used in Eq. the 
fact that 0„(fc) in the two-dimensional system is per- 
pendicular to the plane. Similarly we obtain the edge 
current for the edge along the x direction I^^ which is 
identical to ly. Thus the edge current does not depend 
on the edge direction or the expression for the confining 
potential U (r). Therefore, the magnon moves even along 
the curved edge. Here we should note that in addition to 
the velocity along the edge (i.e. the second term in the 
r.h.s. of Eq.®), there exists a group velocity (the first 
term in the r.h.s. of Eq.©). Because of this group ve- 
locity, a single wavepacket does not go purely along the 
edge. What we have shown is that there is an additional 
velocity along the edge due to Berry curvature, and the 
total magnon edge current is given by Eq. (jSj when all 
the magnons in thermal equilibrium are summed over. 

If the chemical potential p or temperature T varies 
spatially, the thermal Hall effect will occur because the 
magnon edge current no longer cancels between one edge 
and the opposite edge, and a net current will appear. In 
the following we show the details and calculate thermal 
transport coefficients. We focus on the edge current in 
the X direction with small temperature gradient in the y 
direction as an example, and set the coordinate system 
shown in Fig. [TJ Here w is the width of the system and 
a, &i, &2 is defined as U{a) = 0, U{hx) = U{h2) = oo 
and hi < ~w/2 < a < w/2 < b2- The current density is 
obtained by summing up the local current density jx{y) 
and dividing it by the width: 



jx = — dyjxiy) = - 



1 ['''' 1 r 

- I dyjx{y) + — / dyjx 

W Ja W Jt, 



(9) 



Here we defined jx{y) as the following: 



n.k 



3 



This quantity is nonzero when dU{r)/dy 7^ 0, i.e., y ^ 
zLw/2. At these points, 

p{enk + U{r);T{y)) (11) 

^ fp(e„fc + C/(r);r(f)) (y ^ f ), 
\p(e„fe + [/(r);T(-f)) (y ^ -f ). 

Thus Eq. ^ is written as 
jx = -tttY] [ de 

w hv ^ . 

n.k ^"'= 



= i\wll r p{s:ny))n„Am] ■ (12) 

The edge current along the y direction with temperature 
gradient in the a; direction is similarly written as 



3v 



d_ I J_ 



J2 r p{e;T{x))n^,,{k)de \ . (13) 



In the presence of the chemical potential gradient the 
edge current can be written as the same form like 
Eqs. and (IT51) . Therefore, if the spatial variation 
is well gradual, the edge current density is written as: 



1 f°° 
J- Vx — V / p(£)f2„(fe)de 



(14) 



In the same way, the energy current from the edge current 
density is written as 

1 f°° 

JE = V X — ^ / £p(e)f2„(fc)&. (15) 

n,k •^'^"fc 

From Eqs. (jl4l) and (jlSp . we can derive various transverse 
transport coefficients. For instance, in the presence of the 
temperature gradient in the y direction again, the edge 
current and energy current density in the x direction are 
written as 

o)r^™.(i)i:££^(|)o..,.M. 

(16) 

oe)r^™.(i)i:£i^(|)o,.,.,.)... 

^ ^ n,k ""^rik \ / 

(17) 

Here we note that the temperature gradient affects 
these currents through the Bose distribution function as 
drp (e) = (e - m) {dp (e) /de) Tdr (1/T). Similarly we 
obtain these currents in the presence of the gradient of 



the chemical potential in the y direction: 

n.k 



(18) 



dp 



e I — ) nnAk)de. (19) 



Now we define a heat current as j'q = Je — nj, and write 
down the linear response of the magnon current and heat 
current as 



j ^ Lii [-VU - V/i] + L12 
jq = L12 [-VC/ - V/i] + L22 



TV 
TV 



(20) 
(21) 



From Eqs. (IT51) - p^ . the transverse thermal transport 
coefficients Lf^ can be derived as 

^^f = -WmY.^-^^-^^^^M, (22) 

n,fc 

where i,j = 1,2, c,(p„) = j2je{P{e-iAf{^^) = 

Jo" (log dt, q = i+j -2, and p„ ee p(e„fc). For ex- 
ample, cq{p) = p, Clip) = (1 + p)log(l + p)-plog/3, and 

C2(p) = {l + p) {log'-^y - {logpf - 2Li2(-p), where 

Li2(2;) is the polylogarithm function. Finally we derive 
the thermal Hall conductivity in a clean limit by substi- 
tuting Eq. (I22|) to K'^y = Lll/T, 



K 



hv 



.,k 



dkr 



dk„ 



(23) 



Thus the thermal Hall conductivity is expressed as the 
Berry curvature in momentum space, which is sensitive 
to the magnon band structure. Since the Berry curvature 
part is expressed as 



Im 



dH 



dH 



dk-u 



(24) 

Hence in Eq. ([25)1 is enhanced if there is an avoided 
band crossing. 



III. LINEAR RESPONSE THEORY 



Compared with our result in Eq. (j23p . the expression 
for the thermal Hall conductivity obtained in the previ- 
ous works^, 



hVT 



n,k ^ 



H + £nk 



du„ 



(25) 
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is different in some points. As we see later in this sec- 
tion, in the hnear response theory, our resuh (Eq. 
consists of two parts (S')f/ + {M)J^ (defined in Eqs. ([SS]) - 
(p4| ). while the result in Refs. (Eq. ^5^) contains 
only {S)1J . The correction term {M)^J comes from or- 
bital motions of magnons. In the following, we apply the 
linear response theory to the magnon system in analogy 
with the electron systemiir— . In Appendix \^ we briefly 
review the linear response theory under the temperature 
gradient in the electron system, and derive some useful 
expression for the thermal transport coefficients. 

Here we shortly discuss the linear response theory with 
external fieldsii"— . In the presence of the temperature 
gradient, it is convenient to introduce a fictitious gravi- 
tational potential ?/'('") i which exerts a force proportional 
to the energy of the particle^. This is because in order 
to obtain the thermal transport coefficients by the linear 
response theory, we need to take the temperature gradi- 
ent into the Hamiltonian as an external field. However, 
it is not possible to directly incorporate the temperature 
gradient into the linear response theory, since the tem- 
perature gradient is not a dynamical force which exerts 
force to the particles, but a statistical force which affects 
the particles through the distribution function. There- 
fore, to avoid this difficulty, the fictitious potential i/;, 
giving a dynamical force, has been introduced. As we see 
in Appendix [XI the thermal transport coefficients from 
the temperature gradient are derived by calculating the 
coefficient from the gradient of the fictitious potential 
Tp. This is analogous to the situation that the transport 
coefficients from the gradient of the chemical potential 
can be obtained by calculating the coefficients from the 
electric field. 

Now we consider the magnon system. Since the 
magnon has no charge, we cannot use the electric field 
E as an external field. Instead, we again use the gradi- 
ent of the confining potential — VJ7(r), which appeared 
in Eq. The perturbation Hamiltonian is written as 

H' = U{r,) + \{H,j, r, ■ VV^(r)}, where 

the position of the jth magnon, H is the unperturbed 
Hamiltonian, and {A, U} — AB + BA represents the an- 
ticomutator. In equilibrium, the magnon current density 



and energy current density are written as 

(26) 

j 



(27) 



where Vj is the velocity operator of the jth magnon. In 
the presence of the external fields H\ they acquire addi- 
tional terms, 

j(r-) = + ^ |/°Hr),l^r, . VV(r)| , 

(28) 

j 

+ i^E({{^'^^-^VXr-)},j(°Hr)} 

j 

+ {{j(o)(r),r,.VV;(r)},if}), (29) 



where c is the speed of light. Correspondingly, the ther- 
mal transport coefficients consist of two parts: {L)"f — 



iS)tf + {M)IP. Here a,p = x,y, i,j - 1,2, and in 

this section the thermal transport coefficients {L)"/ are 
defined as: 



Je ^ (i)i2 



-Vi7 - TV (J 



-VC/ - TV (^^ 



(i)l2 
(i)22 



TV(i 



7^V( - 



vy; " 

(30) 
(31) 



A deviation of the distribution function from the equi- 
rj is librium state generates {8)°'^^ , which is calculated by the 
Kubo formula; a deviation of the current operator due 
to external fields from the equilibrium state generates 
{M)"f . In the clean limit they are expressed as 



(S^)tf -yj Piv)Tr (jr^J^Siv -H)- jrSi^ H)jf^^ dr,, (32) 
(M^)?f = 0, (M^)?2^ = pivmSiv - HKr^vP - rPv^)]di^, (33) 

(M^)^l - / VPm^Siv - H){t'^v^ - rf^v^drj + ^ f p{ii)TvS{v - H)[v'^ ,v^]d7j. (34) 



r 



Here G'='= is the Green's function G='=(?7) = (r?- i/±ze) ^ distribution function p{ri) = (e'^*''"^) - l) \ ji v, 
with e being the positive infinitesimal, p{ri) is the Bose = ^(^Hv + vH), v is the velocity of magnons, and the 
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label of the superscript "B" means a boson. By using 
Eqs. (|A10p - (|A13p . these thermal transport coefficients for 
the magnon system can be written by the wave function 
of magnons: 



n.k 



H + £nk 



B\al3 
ij 



^^ImX:c,(p„)(^ 



(35) 
(36) 



where we have taken /x = since the magnon num- 
ber is not conserved. We note that q — i + j — 2 and 
(M^)"f = 0. Thus the total thermal transport coeffi- 
cients are written as 



,rB\al3 
ij ' 



(37) 



n,k 



and from this equation we again derive the same thermal 
Hall conductivity as Eq. ([23]). The result in Refs. US 
shown in Eq. (j25p . contains only the contribution from 

{8^)22 ■ Therefore the difference between the results 
of Refs. [^1^ and ours arise from the correction terms 

As we can see from Eqs. (p3| and ((34| . these correction 
terms {M^)"^ are related to the orbital motion of the 
particle, namely, a reduced orbital angular momentum 
(rxv). Equation. (p6|) means that {M^)"^^ are expressed 
as the Berry curvature in momentum space, which is gen- 
erally nonzero. Hence, in this case, the magnon has fi- 
nite orbital angular momentum due to the Berry cur- 
vature. This orbital angular momentum consists of two 
parts: the edge current and the self-rotation motion of 
the wavepacket. The reduced angular momentum for the 
edge current per unit area is derived from Eq. (|8]), 



pdge 



n.k 



(38) 



and that for the self-rotation motion is calculated in anal- 
ogy with the electron systemic as 



self 



2 / dun 



{H - £„fc) 



regarded as a generalized cyclotron motion. However, 
since the magnon has no charge, it feels no Lorentz force 
and cannot have a cyclotron motion in the same sense 
as that of electrons. In this respect, this motion purely 
reflects the magnon band structure. A similar effect can 
be found in various wave phenomena such as electrons^, 
photons'^^, and so on. 

As mentioned earlier, within the semiclassical theory, 
the result for the edge current is derived under the as- 
sumption that the spatial variation of the confining po- 
tential is slow. Nevertheless, as we have seen in this 
section, the linear response theory, which does not need 
assumptions for confining potential, gives the same trans- 
port coefficients as the semiclassical theory. This strongly 
suggests that the edge-current picture carries over to the 
abrupt spatial variation of the confining potential. In 
the quantum Hall systems, this idea is indeed true; at 
the abrupt edges of the quantum Hall system, the elec- 
trons undergo a skipping motion. Namely, near the edge, 
electrons undergo a cyclotron motion and when electrons 
hit the edge they are bounced. As a whole the electrons 
go along the edge with skipping orbitals, which are re- 
garded as the chiral edge current in the quantum Hall 
system. Therefore, we can similarly expect that in the 
ferromagnets with edges, the magnon will undergo a su- 
perposition of a skipping motion along the edge and a 
motion along the group velocity of the magnon. 

The coherence length of the magnons is important in 
the orbital motions and transport of magnons. For the 
validity of the linear response theory developed above, 
it is implicitly assumed that the coherence length of the 
magnons is sufficiently short compared with the system 
size. By this assumption, when we apply temperature 
difference between the two opposite sides of the system, 
we can define a local temperature, and the temperature 
gradient becomes uniform. The linear response theory 
is then justified. Otherwise, when the coherence length 
is as long as the system size, the magnon transport is 
described in the similar way as the Landauer formula, 
and the linear response theory no longer applies. 

Orbital motions of electrons generate magnetic mo- 
ments by their charge. On the other hand, magnons do 
not have charge, but have magnetic moments. There- 
fore we can regard the rotating magnon as a circulating 
spin current. As is similar to the magnetoelectric ef- 
fect in noncoUinear spin structure^S, the rotation of the 
magnon is expected to generate an electric polarization. 
For this effect, the spin-orbit coupling interaction, such 
as Dzyaloshinskii-Moriya (DM) interaction, is necessary. 



It is easy to show that if^" + if^ = 2{M^)^^. This 



result is expected from the Eq. (|33|, i.e., the correction 
terms comes from the orbital angular momentum of the 
magnon. 

Therefore, due to the Berry curvature, the magnon 
generally has a nonvanishing orbital angular momentum 
in equilibrium. This orbital motion of magnon can be 



IV. EXAMPLE 1: LU2V2O7 

In this section we apply our results to the ferromag- 
netic Mott-insulator LU2V2O7 with pyrochlore structure, 
for which the thermal Hall effect has been measured and 
analyzed in Ref. \^ Following Ref. |^, we briefly review 
the magnetic properties of the material. The magne- 
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tization comes from spin- 1/2 V'''''' ions with the DM 
interaction. The ground state is a cohinear ferromag- 
net, because at the ferromagnetic ground state the to- 
tal DM vectors for the six bonds sharing a single site is 
zcrc^. The DM interaction affects the spin- wave disper- 
sion, and the effective spin-wave Hamiltonian is written 



as Hcff = -JSi-Sj+D,j-{Si X Sj)~g^iBH-Y,i Si, 

where denotes the nearest neighbor pairs, J is the 
exchange interaction, D is the DM vector, g is the g- 
factor, /iB is Bohr magneton, and H is the magnetic field 
in the z direction. The temperature is assumed to be 
much lower than the Curie temperature Tc — 70 [K], for 
existence of well-defined Bloch waves of magnons. There 
are four magnon bands, and the lowest band is well sep- 
arated from the other higher bands, with the separation 
much larger than ksT. Actually, the differences of the 
energy between the lowest band a nd other bands near 
k = are written as 62 — Si ~ AJSy/l + f{k) ~ 8JS and 
£3 - £1 = £4 - £1 ~ AJS + 2JS^1 + f{k) ~ 8JS, where 
/(fc) = cos{2kxA)cos{2kyA) + cos{2kyA) cos{2k^A) + 
cos(2fc^A) cos(2fc^A) and 8JS ~ 13.6[meV]. Therefore, 
the contribution from the lowest band is dominant, whose 

Berry curvature is Vli^z ~ "sv^T"^'-^^ + + ^^z) as 
calculated in Ref. 0, with A being a quarter of the lattice 
constant. We can estimate the orbital angular momen- 
tum of the magnon from both the self-rotation motion 
L|°'^ and the edge current L|'^se_ Near k = 0, the lowest- 
band dispersion is quadratic and we can introduce the 
effective mass of the magnon of the lowest band m| , de- 
fined as TO* = h'^{d^enk/dk'^)~^ . The orbital angular 
momentum of the self-rotation motion is analytically cal- 
culated from Eq. (1551) : 



JSml D 1 



WJSml 
W 

knT 



Im^p(£ifc)(|^ 



HA J 247r2 
JSml D 1 



5/2 



dka 

3/2 



dui 
dkp 



JS J JO 

5/2 



dx 



hA J 327r3/2 \js) 



Lis 



(40) 



We obtain ~ -0.009ft and if sc ^ mllf'^" ~ 

-|-0.008/i per unit cell. The thermal Hall conductivity 
is also calculated, by assuming that the contribution 
of the lowest band dominates. Figure [5] show the result 
of the thermal Hall conductivity which is calculated from 
{S)'^l^ + {M)"^^ (solid curve) and {S)ff (broken curve). 
They correspond to our results and the previous results 
in Ref. 6, respectively. Our result (solid curve in Fig. [2]) 
roughly agrees with the experimental data in Ref. 



EXAMPLE 2: 



MAGNETOSTATIC SPIN 
WAVE 



In the following, we apply our theory to the magneto- 
static spin waves in a ferromagnet. In the magnetostatic 
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FIG. 2: (Color online) Dependence of the thermal Hall con- 
ductivity on a magnetic field. The red (broken) curve denotes 
the result which is calculated from only {S)"^ calculated in 
Ref. y; the green (solid) curve denotes the result which is 
calculated fi-om (S')f/ + (M)^"/. 
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e / 










YIG film 

FIG. 3: Geometry of the coordinate axes. The magnetization 
M precesses around Mq. 



spin wave, the wavelength is sufficiently long and the 
exchange coupling between spins is negligible. The mag- 
netic anisotropy comes from the demagnetizing field de- 
termined by the sample shape. This magnetic anisotropy 
due to the demagnetizing field plays the similar role as 
the spin-orbit coupling in electronic systems. It then in- 
duces the Berry curvature, and the Hall effect of spin 
waves appears. 

Let us consider the yttrium-iron-garnet (YIG) film 
which is magnetized by an external magnetic field. We 
introduce two coordinate systems xyz and ^77^, shown 
in Fig. |3l The film is taken to be infinite in the 77- and 
(■-direction, and perpendicular to the ^ direction. C axis 
is chosen to be along the magnon wave vector k. The z 
direction is parallel to the saturation magnetization Mq 
and the internal static magnetic field Hq. We assume 
that the spin wave mode has a form of the plane wave: 
"^(C: Cj t) = "T'(C) exp(i(fc( — wt)), where w is a frequency 
of the spin wave . The equation of motion of the magne- 
tization is written as the following integral equation^i: 

rL/2 

ujHm{0-ujM / dme,^')M^')=waj,m(0. (41) 

J-L/2 

Here we use the SI units, L is the thickness of the film, 
(Ty = the PauU matrix, m{£) = (^^[5]) is 

the vector Fourier amplitude of the spin wave which is 
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perpendicular to Mq, luh = jHq, ojm = iMq, and 7 
is the gyromagnetic ratio. G'(^,^') is the 2x2 complex 
matrix of the Green's function: 



(42) 



^{^jt. ) ~ \ r' r / ' 

Gxx = (Gp - S{C - C')) sin^ ^ - iGg sin 29 cos 

- Gpcos^e'cosV, (43) 

Gxy Gyx 

= —iGq sm9siinp — Gp cos0sin</5cos(p, (44) 



Gyy = —Gp sin (fi, 
where 



Gp = -exp(-fc|e-e'|), 
Gq = Gpsign(e - a, 



(45) 

(46) 
(47) 



and 6, ip are the spherical coordinates of Mq in the ^77^- 
space (see Fig. [3]). We note that we adopt the defini- 
tions of 6 and ip used in Ref. 21, and they are differ- 
ent from the standard definition of the spherical coordi- 
nates. The integral equation (pij) is equivalent to the lin- 
earized Landau-Lifshitz equation dM/dt = —j{MxH), 
Maxwell equation in the magnetostatic limit V x Jf = 0, 
W ■ B — 0, and the usual boundary conditions for H 
and B. Since the equation (PT|) is a generalized eigen 
value problem, we have to modify the prescription of our 
theory of the Berry curvature. Similar to the previous 
worls^^, the Berry curvature is defined as 



^n,7 (fc) 



drrin., 



dka 



dnXr. 



(48) 



Since all the terms in Eq. (|49|) are real, the eigen vector 
m! is also real. Correspondingly, the Berry curvature 



Eq. (j48|) becomes 



^^n,7 (fc) 



dm' 



dka 



U-^ayU 



dm' 



dkp 



am' 



dkr, 



(53) 
(54) 



Because m' is real and there is no imaginary part, this 
Berry curvature vanishes. Thus when Mq is in the film, 
we cannot expect either an orbital rotational motion of 
spin wave packet or the thermal Hall effect due to the 
Berry curvature effect. In other words, in the magneto- 
static backward volume wave (MSBVW) and the mag- 
netostatic surface wave (MSSW), the effects of the Berry 
curvature do not appear. 

On the other hand, the Berry curvature is finite if 
the saturation magnetization is perpendicular to the film 
{9 = 0), i.e., in the magnetostatic forward volume wave 
(MSFVW). In the following, we demonstrate the calcu- 
lation of the Berry curvature for MSFVW. We note that 
^ coincides with z direction when 9 = 0. The solution of 
the integral equation (|41l) of the n-th band for 61 = is 
written as^ 



cos 



(55) 



where Sap-y is the antisymmetric tensor, n is the band 
index of the spin wave mode, and the bra-ket product 
means an usual inner product of vectors and integral over 
z. 

In some cases, Eq. (PSI) becomes zero because of the 
symmetry of the system. This occurs when the saturation 
magnetization Mq is in the film {9 = 7r/2). When 9 = 
7r/2, we can show r2„ .y(fc) = explicitly by performing 



a gauge transformation m' = U ^m = 



1 
i 



m. Then 



Eq. (1411) becomes a generalized eigen value problem with 
real coefficients: 



cunm'iO-^M I de'G'(^,r)m'(C') 

-L/2 



-uj(Txm'{£,), 
(49) 



where G'(C,C') is 



G'{^,0 = U-'G{i,e)U 

'yx ^yy/ 

G'xx = Gp-5{^-a, G' =-GQsin^, 



G 



'f 

XX 



gsm(p. 



GL 



-Gp sin Lp. 



(50) 

(51) 
(52) 



where k = ujm(^h / {(^h - ^^1), v = ojMUJn/{^H P ■ 

— 1 — K > 0, LOn is the n-th band energy for n = 0, 1, 2, . . 
which is determined by 



■y^tan 



= 1, 



(56) 



and A'^ is a normalization factor which is determined by 



{rrinMl <^y 



(57) 



To obtain Eq. (|55p. we have rewritten the solution in 
Ref. [2^ in the polar coordinate into the form of the plane 
wave. The dispersion determined by (j56l) is shown in 
Fig. Ufa) for Hq/Mq = 1.0. We use the normalization 
Eq. ([57| . because m\ f^cjymn.k is proportional to the en- 
ergy density for the magnon^"*. 

Substituting this solution to Eq. pS)) . we can can 
calculate the Berry curvature for the n-th MSFVW 
mode. For simplicity of the notation, we set _F„(fc, z) = 
\/N cos [y/pkz + ^) . Then the Berry curvature is writ- 
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ten from Eq. (|48 



ri„,.(fc)/2 



= Re 



dmi 



dky 



dky 



drrvi 



dkx 



J-L/2 



dn dv 

k K—- + 



dk 



L/2 



L/2 



dk dk 

Using the normalization condition of Eq. ([5 



Derivative of Eq. (|5. 



2Kvk^ / dzF^ = 1. 

J -L/2 

in terms of k leads to 



dzFl (58) 
one obtain 
(59) 



dzFn - 

-L/2 



dk 



1 dK 
2^dk 



Thus Eq. (|58p is rewritten as 



1 dv 
2^dk 



1 dn 
K dk 



1 dv 
V dk 



L/2 



L/2 



Since k and v satisfy the following relation 

1 9k 1 dv 1 dujn 
K dk V dk uj„ dk ' 

the Berry curvature is derived as 



F^^dz. 
(60) 

(61) 
(62) 
(63) 



Figure |31[b)-(d) shows the numerical results of Eq. 
for various magnitude of the magnetic field. It is surpris- 
ing that the Berry curvature for any MSFVW mode is 
always positive, because ujh < w„ and the group velocity 
duj/dk is positive. In the vicinity of /c = 0, we can cal- 
culate asymptotic forms of the Berry curvature. When 
fc ~ 0, is close to loh- If we set w„ = ujh + Acj„, p 
can be written as p ~ Wj\//2Acj„ ^ 1 since Aaj„ is small 
near fc = 0. Using an approximation tana; ~ x (a; <C 1) , 
we find 



ujm f kL 



{n = 0) 
{n > 0) 



(64) 



Therefore, the Berry curvature near fc = can be ob- 
tained from Eq. and (IMl) : 




{n = 0) 



(feL)2 (n > 0) 



(65) 




10 20 30 40 



10 20 30 40 



FIG. 4: (Color online) (a) Dispersion relation for the MSFVW 
mode with n = 0, 1, . . . , 5, and the Berry curvature for the 
MSFVW mode for (b) Ho/Mo = 0.1, (c) Ho/Mo = 1.0, and 
(d) Ho/Mo = 2.0. 



It is easy to see that f7„_2(fe = 0) = for n > mode, 

and that for n = mode flo,z{k) enhances up to (^jj^ 
but does not diverge at /c = 0. 



VI. CONCLUSIONS 

In summary, we found that magnon wavepacket has 
two types of orbital motions due to the Berry curvature 
in momentum space: the magnon edge current and the 
self-rotation motion. The magnon edge current causes 
the thermal Hall effect of magnon, and the self-rotation 
motion of magnon without Lorentz force is expected to 
accompany an electric polarization. We showed that our 
theory is applied to not only the exchange spin wave 
(quantum-mechanical magnon) e.g. in LU2V2O7, but also 
the classical magnetostatic waves e.g. in YIG. In both 
cases, the Berry curvature is enhanced near the band 
crossings, where the magnon frequency in a focused band 
is close to those of other bands. We expect to control the 
Berry curvature by designing magnonic crystals2^. 
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Appendix A 

Here we briefly review the linear response theory for 
the electron system with a temperature gradient, devel- 
oped in Refs. IllMlSl In equilibrium, the electric current 
operator j^'^\r) and the energy current operator j'^\r) 
are given by 



j^'\r) = --J2{vj,S{r~r,)}, (Al) 

j 

3^EHr)^-l-{H,j('\r)}, (A2) 

where rj denotes the position of the jth electron, — e 
(e > 0) is the electron charge, Vj is the velocity operator 
of the jth electron, H is the unperturbed Hamiltonian 
of the system, and {^4, B} = AB + BA denotes the an- 
ticomutator. We note that rj is a quantum mechanical 
operator, while r is a c-number. 

Under the electric potential (j>{r) and the fictitious 
gravitational potential ip{r), the Hamiltonian is written 

as Htot =H + eY.jr,-V<j,+ \ [h, ^ r, ■ Vi>{r)], 

where H is an unperturbed Hamiltonian and c is the 
speed of light. Subsequently, the electric current opera- 
tor j{r) and the energy current operator jE{f) deviate 
from the equilibrium state. They are written as: 



:jW(r) + i|j(°)(r),l^r,.V^(r)|, 

(A3) 



AO) 

JE 



3 

{{/°'(r),r,-VV'(r)},i/}). 



(A4) 



The linear response for the electric current and energy 
current is written as 



Je = (i^)l2 



(i^)l2 



TV 



TV 













(A5) 












(A6) 



where "F" is a label which means a fermion, E is an 
electric field, /i is the chemical potential, and {L^)ij is 
the transport coefficients — 1,2). The measurable 
current densities J and Je are obtained by taking av- 
erage over the volume of the sample and the quantum- 
mechanical and thermodynamic averages of the current 
operators j{r) and jE{r), respectively. Due to the devi- 
ations j'-^-*(r) and J^^''(t), the thermal transport coeffi- 
cients {L^)'j'f (a,/? = x,y) consist of two parts, {S^)°'f 
and [Antf: 



Here is the Green's function G^{r]) — [r] — H ± 
ie)~^ which is introduced in Eq. (jA7p via S{r] — H) = 
— (G^ — G^)/27ri, f{-q) is the Fermi distribution function 



f{ri) = [e^i^-t-) + 1) \ ji^ -ev, 32 ^ \{Hv + vH), 
and V is the velocity of electrons. {S-^)'^f is calculated 
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from the current operators in equilibrium state, j^°\r) 
and j''^\r), with the deviation of the distribution func- 
tion from the equihbrium state; {M^)'^^ is calculated 
from the deviation of the current operators, j'^^^ir) and 
''('') J with the equilibrium distribution function. Actu- 



ally, {S^)"f is the Kubo formula, and {M^')"f represent 



F\aP 



correction terms. The total thermal transport coefficients 
arc their sums: = {8^)1^ + {MP)f . 

From these results, we can derive some useful equations 
for later calculations. First, we can write down (S^)"^^ in 



terms of the Berry phase. For example, (S*^)"^ is written 
as: 



12 



e 

W 



Im^/(, 



dka 



dkp 
(AlO) 



Second, because the Fermi distribution function /(r/) be- 
comes the step function Q{ii — r]) in the zero temperature 



limit, {L^)i2 and {S^ is written as 



F\aP 



(All) 



and 



H + £„fc 



X 



dUn 

dka 



dkj3 



^ • (A12) 



^^"'3 and 



Thus the relation (L^ = (S" )i2 + {M^ )°2 
Eq. (|A8[) in the zero temperature limit lead to the fol- 
lowing useful formula: 



Tr[(5(^-i?)(r"i;^ -r'^-y")] 



)]dv 



2V_d_ 
—e dfi 
2 d 



(M 



F\al3 
12 



T-J-O 



X Im 



dUn 

dk„ 



[H + Enk - 2Ai) 



dkp 



«P ■ (A13) 



We note that this equation does not depend whether the 
particles are fermion or boson. Therefore we can apply 
this equation to the magnon system as well. 
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